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Abstract 

We study the properties of the Green's functions of the fermions in charged Gauss-Bonnet black 

hole. What we want to do is to investigate how the presence of Gauss-Bonnet coupling constant 

^— I a affects the dispersion relation, which is a characteristic of Fermi or non-Fermi liquid, as well as 

O what properties such a system has, for instance, the Particle-hole (a)symmetry. One important 

{>Q result of this research is that we find for q = 1, the behavior of this system is different from that of 

the Landau Fermi liquid and so the system can be candidates for holographic dual of generalized 

non-Fermi liquids. More importantly, the behavior of this system increasingly similar to that of 

r"| the Landau Fermi liquid when a is approaching its lower bound. Also we find that this system 

» 

possesses the Particle-hole asymmetry when (? 7^ 0, another important characteristic of this system. 

(D 

4^ In addition, we also investigate briefly the cases of the charge dependence. 
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I. INTRODUCTION 



As is known to all, the success of the single electron picture of metals rests on Landau's 
Fermi-liquid theory. The metal in such theory is treated as a gas of Fermi particles whose 
interactions are weak and not as significant as that of the original electrons. The reason is 
that particles of this model are not the original electrons but the electron-like quasi-particles 
that emerge from the interacting gas of electrons. However, recently some new materials, 
including the cuprate superconductors and other oxides, seem to lie outside this framework. 
For these new materials, we refer them to the non-Fermi liquid metals. They show lots 
of new physical properties which can not be understood in terms of weakly interacting 
electron-like objects. For the non-Fermi liquid, a sharp Fermi surface still exists. But quasi- 
particle picture breaks down generically. Although there have been many phenomenological 
models to describe the non-Fermi liquid, a general theoretical framework characterizing 
non-Fermi liquid metals remains a suspense. Therefore, it is necessary to develop a basic 
principle for non- Fermi liquid. Maybe the AdS / CFT correspondence can provide us with a 
possible clue to yield the basic principle of non- Fermi liquid. Indeed, by applying AdS/CFT 
correspondence P-[3], some breakthroughs have been achieved, that is, we find some new 
classes of non- Fermi liquids |4j. Varieties of holographic fermions models and their extensions 
are being explored. These can be seen from Refs. [5]-[12]. 

However, before string theory is fully understood, it is necessary to consider the higher 
curvature (or derivative) interactions in an effective gravitational theory. From the point 
of view of the AdS/CFT, the higher curvature interactions on the gravity side correspond 
to finite coupling corrections on the gauge theory side, thus broadening the class of field 
theories one can holographically study. The main motivation of considering these corrections 
derives from the fact that string theory contains higher curvature corrections arising from 
stringy effects. A simple and useful model with regard these corrections is Gauss-Bonnet 
(GB) gravity, which contains only the curvature-squared interaction. Several models of 
holographic superconductor have been exploited in this setting [T314T9] . Recently, a new 
higher derivative theory of gravity (quasi-topological gravity) is constructed, which contains 
not only the curvature-squared interaction but also a curvature-cubed interaction [201 |2T] . 
The corresponding holographic superconductor models have also been discussed in |22l |23] . 
Another interesting extension comes from the coupling between the Maxwell field and the 
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bulk Weyl tensor [2ll |25]. Also, in Refs. [251 126], they find that the Weyl corrections can 
describe the dual field theories with both the weak-coupling and the strong-coupling, which 
is a very interesting and significant holographic study on the dual field theories. Therefore, 
it is deserving to further exploit the effect of higher curvature interactions on fermions by 
using AdS/CFT correspondence step by step. In this paper, we consider the Gauss-Bonnet 
term as the first step to introduce the stringy correction into the gravitational action. And 
we will mainly focus on how the GB coupling constant affects the spectral function of the 
fermions. 

Our paper is organized as follows. In section II, we briefly introduce the charged Gauss- 
Bonnet black hole. Following Ref. j^, we obtain the Dirac equation of the probe fermions in 
the charged Gauss-Bonnet black hole in section III. The numerical results are presented in 
section IV where we will mainly focus on the dispersion and the Particle-hole asymmetry in 
particular. Conclusions and discussions follow in section VI. Finally, in appendix A, we give 
a brief summary on the analytic treatment in the low frequency limit developed in Ref. [H] . 



II. CHARGED BLACK HOLES IN GAUSS-BONNET GRAVITY 

In this section, a brief review of black holes in Gauss-Bonnet gravity is given as follows. 
Conventionally, we begin with the following action coupled a vector fleld Aa ^: 



S 



R 



abed p 

^abed 



ab 



R^) 



pabp 



ab 



(1) 



where Rabedi Rab and R are the Riemann curvature tensor, Ricci tensor, and the Ricci scalar, 
respectively. L is the AdS radius and for convenience, we will set L = 1. a is Gauss-Bonnet 
coupling constant and the constraints on it will also be discussed in the following parts of 
this paper. And as is commonly known. Fab = daAb — dbAa- By applying the principle of 
variation to ([T]), we can easily obtain the equations of motion: 
Maxwells equations 



ab 



0, 



(2) 



^ Here we have set the the five dimensional gravitational constant Kg = 1/2 and the effective dimensionless 
gauge coupling gp = 2. 
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and Einsteins equations 



Rab - ^Rgab + Qgab - a 



Hab — -^Hgab 



abi 



where 



and 



TJ rycde p o D DC o p rycd I 75 D 

J^ab — J^a J^bcde — -^^acJ^b ~ ^J^acbdii- + rCKab- 



H = m. 



(3) 



(4) 



(5) 



Under a certain circumstance in which we take the following metric with planar symmetry 



ds^ = -f{r)dt^ + jj-^ + r^(dx^ + dy"^ + dz^), 



and the ansatz of the gauge fields as 



A„ = (A(r), 0,0,0), 



the equations of motion ^ and ^ reduce to 



+ -^t = 0, 
r 



+ -/ - 4r + — ^ = 0, 
r D 



(6) 



(7) 



(8) 



(9) 



where the prime represents derivative with respect to r. The solutions of the above equations 
are [271 128] 

^2 



fir) 



2a 



1 - W 1 - 4a 1 



+ 



3r4 



and 



At = fill 



,2 y ' 



(10) 



(11) 



where r+ is the horizon radius determined by f{r^) = 0^, and fi can be identified with the 
chemical potential of the dual field theory. In the Einstein limit a — )■ the formula above 
reduces to the Reissner-Nordstrom AdS black hole 

,2 



/(r)Uo = r^ 



3r4 



(12) 



In the following, for convenience, we will set ?'+ = 1. 
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By using the standard approach of euchdean continuation near the black hole horizon, 
the Hawking temperature of the black hole in Gauss-Bonnet gravity is 

(13) 



Ait 71 \ 6 ^ 

which is also the temperature of the conformal field theory on the boundary of the AdS 
spacetime. It is independent of the Gauss-Bonnet coupling constant a. 

In the following, we will discuss the behaviors of places verging on the boundary (r — )■ oo) 
and the horizon (r — )■ 1), respectively. Near the boundary of the bulk, the redshift factor 
/(r) becomes 

2 

/oo = lim /(r) = ^ [1 - VT^] = r^Ca, (14) 



where we denote = (1 — \/l — 4a)/2a. Then, the geometry of the black hole on the 
boundary can be reexpressed as 

ds"^ = -r^Cadt'^ + + r^{dx^ + dy"^ + dz'^). (15) 

It is also pure AdS^, but depends on the GB coupling constant a. In order to have 
a well-defined anti-de Sitter vacuum for the gravity theory, a certain condition a < 1/4 
should be required, whose upper bound a = 1/4 is known as the Chern-Simons limit. 
Furthermore, considering the causality of dual field theory on the boundary, there exists a 
stronger constraint on the GB coupling in five dimensions [291 - 138] 

7 9 

- — < a < — . (16) 
36 - - 100 ^ ^ 

As for the condition of the horizon, we in reality do a similar job, but for this part, 
we consider a special case, that is, the zero-temperature which means fi = Then, the 
redshift factor /(r) becomes 



12a 8a 

4a H 



(17) 



Obviously, when r — 1, /(r) ^ 12(r — 1)^, which is also independent of the Gauss-Bonnet 
coupling constant a. In light of this, near the horizon, the geometry is no longer the pure 
AdS^, but the AdS2 x with the curvature radius of AdS2, L2 = -4^L. 
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III. DIRAC EQUATION 



Now, we consider probe fermions in the charged Gauss-Bonnet black hole. We have the 
bulk fermion action as following [1] 



where is the covariant derivative given by^'^ 



(18) 



^iMu)a is the spin connection 1-forms given by 



(19) 



(20) 



and 



2^ ' 



(21) 



where (e^)*^ form a set of orthogonal normal vector bases. The Dirac equation derived from 
the action Sd is expressed as 



We should choose the following orthogonal normal vector bases 

, d 

(e^)" = v^(i^)", f^ = t,x, y, z, r. 



(22) 



(23) 



Using (20), one can calculate the non- vanishing components of spin connections as follows 



i(^ir)a = -{(^ri)a = -\f(rdr{y/(hi){dxla, i = X,y,Z. 



(24) 



^ Throughout the paper, we will use the conventions of |39j . a, b are the usual spacetime abstract index 
and fj,, V are the tangent-space index. 



Note that the Dirac action (18 1, only through the effective chemical potential = iiq, that is to say, 
through the combination of gpl, depends on q. For convenience, the gp has been set as 2 above and the q 
is treated as a free parameter, but we should note that only the product of them is the relevant quantity. 
For more discussions, please refer to Ref.[l]. At the same time, we also want to remind readers to pay 
attention to the fact that in Ref . [H HT] . they set gp = 1, which is different from our conventions. 
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Following Ref.[^, we can make a transformation ( = {—gg^^) to remove the spin 
connection in Dirac equation. Then, the equation turns out to be 

y/g^rdrT + V7*r*(5t - iqAt)^ + ^/g^^^^^ T - m^T = 0. (25) 

Next, we will work in Fourier space where we expand J-" as J-" = ^e-«'^*+«'=i^\ Then, the 
Dirac equation can be rewritten as 



^/^VdrF - i{uj + q) ^/g^^rF + iky/^V'F - nn^F = 0. (26) 

where due to rotational symmetry in x — y — z directions, we set ky = kz = and kx = k 
without losing generality. We will choose the following basis for our gamma matrices as 
ingnilil] 

Splitting the 4-component spinors into two 2-component spinors F = (^1,^2)"^, we have 
a new version of the Dirac equation as 

VF^dr \ ^\+ rrii-a^ ® | = v^(^ + ^A)^^^^ ® f T k^^a^ ® I > (28) 

\F2j \F2j \F2j \F2j 

a decouple equation between Fi and F2. After achieving this equation, we will discuss it in a 
special case, that is, near the boundary (r — )■ 00) in order to acquire some information about 
the Green function near the boundary. From the forgoing information, we know that in this 
special case, the geometry is an an asymptotic AdS^, grr ~ and gu ^ r"^, i = t,x,y,z. 



Under such a condition, Eq. (28) becomes 



{rdr + m^a^) ® f = , (29) 

\F2j 

another decouple equation between Fi and F2 whose solution can be expressed as 



F/^"a„r"^^ (°| +6«r-'"c I ' I a = 1,2 . (30) 



7 











If 1 




and aa 













are related by 





^1^ 






ba 1 




= Saa 1 













(31) 



then the boundary spinor Green functions G is given by 

G = -i57° . 



(32) 



where 7° is the gamma matrices of the boundary theory and 7° = ia^. In order to find the 
matrix S, we will make such a decomposition F± = ^{1 ± T'^)F according to eigenvalues of 
r''. Then 



B2 



A, 
A2 



(33) 



with 



Fa 



A, 



(34) 



Now, by using (30), (31) and (32), we can express the boundary Green functions as 
following 



G{uj,k) = lim r2'"cG(r,w,fc), 

r—^oo 



(35) 



here we have defined the following matrices: 

G{r, u, k) = 

and 



Gil 
G22 



(36) 



Aa 

Gaa{r,u,k) = ,a = l,2. 



(37) 



Moreover, under the same decomposition, the Dirac equation (|28|) can be rewritten as 

kVg^ I I , (38) 



± m^) I = ±{u + qAt)^fJ^ \ 

i3i/ \a. 



Bi 
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A2 
Bo 



, ,E2\ , \B2 

A2] \A2 



(39) 



Using the method developed in 



one can package the Dirac equation (38) and (39) 



into the evolution equation of Gir^ w, k\ which will be more convenient to impose the initial 
conditions at the horizon and read off the boundary Green functions, 



+ 2m^)G' = G [^J~(f{u} + qAt) + ky^a^^ G + {^/g^^iu + qAt) - k^/^a^) .(40) 



The boundary condition of the matrix G(r, cu, k) in this new equation 



G„„(r, k) ''=^ i. 



(41) 



can be derived from the requirement that the solutions of Eqs. (38) and (39) at the horizon, 
r — 7- 1, be in-falling. 

However, we also note that for T = and w = 0, the boundary condition has to be 
modified as follows 



Gaa{r,U = 0, k) 



r— !-l 



(42) 



IV. PROPERTIES OF SPECTRAL FUNCTIONS 



A. General behavior 



As Ref. [4J, we also have some symmetry properties of the Green function by direct in- 



spection of the equation (40). We itemize them as follows: 



(1) G22{uJ, k) = Gn(u;, -k); (2) ^22(0;, k- -q) = G\,{-u, k; q); 
For the case = 0, 



(3) G22itO,k) 



Gii(o^,fc)' 



(4) G22iuj, k = 0) = Gn{co, k = 0) = i. 



When the background geometry is pure AdS^ (Eq. (15)), for massless bulk fermion, the 



Dirac equation (40) can be explicitly expressed as 



+ 



G, 



1/2 



W + yUg k 
Gn 



r 



1/2 • 



(43) 



22 



W + /ig k 
Gn 



G22 + 



Gn 



(44) 
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where we denote the Green function in pure AdS^ background as Q. The solution of the 
above equations can be easily obtained as ^ 



/II 



A 



1 /2 

{U + flq)/Ca - k/Ca +16 



122 



\\ (a; + /i,)/C„-A;/Cy' + 2e^ 



(45) 



{uj + ii,)/Co^ + k/Cl!^ + ie 

where e — ?■ 0. For g = 0, it is clear that the spectral function has a Particle-hole symmetry 
(symmetry under (0;,^) {—uj,—k)). However, when q ^ 0, the Particle-hole symmetry 
is broken (Particle-hole asymmetry). In addition, the spectral function ImQ has also an 
edge-singularity along u = ±k and vanishs in the region u G {—k, k). 

Now, we turn to the charged Gauss-Bonnet black hole background (Eqs. ^ and (10)). 



In this case, the Dirac equation (|40|) becomes 



^ 'u + qAt) + k- 



C,^22 



fl/2 



{u + qAt 



[00 



qAt) - k- 



^22 + 



j^{u + qAt) + k-. 



(46) 



(47) 



We can solve the above equations numerically with the boundary conditions (41) to 



investigate the properties of the spectral function. In this paper, we will only focus on the 
massless fermion (m^ = 0) and extremal charged GB black hole (zero temperature limit). 
The dependence of the spectral function on mass and temperature will be discussed in the 
future works. 

Here we will do several quick checks on the consistency of our numerics. In FIG. [TJ we 
show the spectral function /mG'22 at A; = 2.2 < fig (left plot) and k = 4.0 > fig (right plot) 
for = and q = 1 (/iq = a/6). Firstly, the divergence in the vacuum turns out to be 
a peak of finite size at w + /i^ ~ ±k. In addition, independent of the parameter a, for a 
fixed large k ^ ^g, /mG'22 is roughly zero in the region u + fig E (— /c, k) and asymptote 
to 1 as |c(;| — )■ 00, which recovers the behavior in the vacuum. It is consistent with the 



Green function Q in pure AdS^ background (Eq. (45)). The height and width of the peak 
vary with the GB coupling constant a. When a increases, the peak becomes sharper and 
narrower. We also note that for k = 2.2 < fig (left plot), the deviation from the vacuum 
behavior becomes significant. For some more concrete investigations, we will discuss some 
specific properties of the spectral function in the subsequent subsection. 



^ We can also refer to Refs.f 
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FIG. 1: Spectral function ImG22 at k = 2.2 < fig (left plot) and k = 4.0 > /ig (right plot) with 
different Gauss-Bonnet coupling constant a (a = 0.09 for red line and a = —0.19 for blue line) for 
= and q = 1 {nq = VQ). 

B. Fermi surface and the dispersion relation 

In this subsection, we will focus on the dispersion relation between small k = k — kp and 
u, a characteristic of Fermi or non-Fermi liquid. In order to achieve this goal, firstly, we have 
to find the Fermi surface. As is known to all, the fermion is created near the Fermi surface 
and so it should have a long lifetime. Therefore, when the energy equals the Fermi energy 
and the momentum equals the Fermi momentum {kp), the spectral function of this system 
should have a sharp quasi-particle peak. In this paper, we will adopt the conventions |ll[5], 
where energy equal to the Fermi energy corresponds to the frequency vanishing. Therefore, 



in order to find the Fermi surface, we can solve the equations (46) and (47) numerically 



with the boundary conditions (41). For definiteness, we firstly focus on the cases of g = 1 
(with Hg = \/6) in this subsection and consider the cases of charge dependence in the next 
subsection. For a = —0.19, we obtain a sharp quasi-particle-like peak near kp = 2.071564 
(Fig|2] and Figjs]). Similarly, we can also get different fermi momentums corresponding to 
different a, such as kp ^ 1.7821 for a = 0.09, kp ^ 1.8770 for a = 0.01, kp ^ 1.8879 
for a = 0.0001 and kp 1.8880 for a = (Figjsj). We note that the fermi momentum is 
a-dependent. 

Now, we can move on to investigate the behavior of /mG'22 in the region of small k = k—kp 
and u. By fitting the data, we find that there exists a dispersion relation between A; (A; — t- 0_) 
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FIG. 2: 3d plot of I'mG22{^, k) for a = —0.19, m^^ = and q = 1 {jiq = ^f^)■ A sharp quasi- 
particle-like peak occurs near kp ~ 2.07, indicating a Fermi surface. 



and uj{k) (Fig|4|, i.e., 



Co{k) ~ P 



(48) 



where 5 ^ 1.13 for a = -0.19, 6 ^ 1.35 for a = and 6 ^ 1.55 for a = 0.09. 

In addition, we can also find that the scaling behavior of the height of /mG'22 at the 
maximum as follows: 



(49) 



for all a {FigQ. 

Here come some comments on the above two scaling behaviors. They are different from 
the landau Fermi liquid, which has exponents 6 = P = 1. Therefore, the system can be 
candidates for holographic dual of generalized non-Fermi liquids. However, we also note that 
when the Gauss-Bonnet coupling constant a is approaching the lower bound, the parameter 
6 decreases, indicating the behavior of this system is more similar to that of the landau 
Fermi liquid. It seems that the GB term can also describe the dual field theories with both 
the weak-coupling and the strong-coupling. 

In fact, when the Fermi momentum kp has been determined numerically, the scaling 



exponents 5 in the dispersion relation (48) can also be computed by the analytical method 
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lm(G2j) 




FIG. 3: The plot of ImG22{k) for a; = —0.000001 with different a (blue line for a = —0.19, green 
for a = and red for a = 0.09). Here, q = l {fiq = \/6). 

developing in the Ref.jH]^. We should compare the numerical result with that obtained by 
the analytical method. The result summarized in Table |Tj From Table |T[ we can see that the 
numerical result is agree well with the analytical that. 



a 


-0.19 





0.09 


5 (numerical result) 


1.12612 


1.35345 


1.55103 


6 (analytical result) 


1.14423 


1.38474 


1.59727 



TABLE I: The scaling exponents 5 for different GB parameter a numerically and analytically. 



C. Charge dependence 

As observed in Ref.jU HI], the Fermi momentum kp increases as we amplify the charge 
q. These features are still preserved in the charged Gauss-Bonnet black hole for fixed GB 
parameter a. We list the values of the Fermi momentum kp for a few other values of charge 
q for the primary Fermi surface in Table [iTf . 

After the Fermi momentum kp is determined numerically, we can compute the scaling 



exponent 6 by using Eqs. (Al) and (A2). The results are showed in Table III For compar 



ison, we also present the numerical fitting for q = 1.2 and a = —0.19, 0.09 in Figjsf. From 



^ For details, please refer to the Ref. "51]. We also give a brief summary in Appendix A. 

^ Since the peak of 6*22 becomes sharper as the charge q increases, we can pin down the Fermi surface to 

more digit for larger q. 
® Another scaling exponent /3 is still approximate to 1 independent of a and q. 
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FIG. 4: Left plots: The dispersion relation between k and cu, where k = k — kp and uj{k) is the 
location of the maximum of the quasi-particle-like peak. Right plots: The relation between the 
height of ImG22 at the maximum and k. From top to bottom, a = —0.19, and 0.09, respectively. 
We fit the data as follows: w « -0.855399(-A;)i-i26i2 and ImG22{i^M ^ 75.7802(-fc)-0-9S9264 
for a = -0.19, -1.59304(-A:)1-35345 ^nd ImG22{(^,k) ^ 31.9054(-fc)-°-9^"3'^ for a = and 
Cj ^ -2.21699(-A;)i-55i°3 and ImG22{^,~k) ^ 19.2059(-A;)-i-°°533 ^j. ^ ^ o.09. Here, q = 1 
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Table HI and FigjS} we can see that for fixed GB parameter a, with the increase of charge 
g, the scahng exponent 6 decreases rapidly and will asymptote to 1 for larger values of q. It 
is in agreement with that found in Ref. [1] . When the charge q is increased to certain values 
(for example, q = 1.2), we find that the scaling exponent 6 decreases with the decrease of 
the values of GB parameter a and asymptote to 1 for smaller a. While for larger charge q 
(for example, q = 1.5), the scaling exponent 5 ~ 1 independent of the a. For furthermore 
exploration on how the charge q and GB parameter a affect together the scaling exponent 
5, we leave it for future work. 





q = 0.5 


q = 1.2 


q = 1.5 


a = -0.19 


0.87 


2.581414 


3.3570706 


a = 


0.81 


2.349913 


3.0572678 


a = 0.09 


0.77 


2.212214 


2.8735578 



TABLE II: The Fermi momentum kp for different charge q and GB parameter a. 





q = 0.5 


q = 1.2 


q = 1.5 


a = -0.19 


3.41726 


1 


1 


a = 


4.38389 


1.06558 


1 


a = 0.09 


5.68267 


1.22051 


1 



TABLE III: The scaling exponent 6 for different charge q and GB parameter a (analytical result). 



D. Particle-Hole (A) symmetry 

As mentioned in subsection A, when the background geometry is pure AdS^ and q = 0, the 
spectral function has a Particle-hole symmetry, but when q ^ 0, the Particle-hole asymmetry 
presents. The similar case does occur in the background of charged Gauss-Bonnet black hole 
(Figjoj). In Fig|6| the last two panels below show that the spectral function behavior is not 
symmetrical about the Fermi point (Particle-hole asymmetry) for a = —0.19 and q = 1 
or q = —1. Moreover, we also find that when q is restored to 0, the system regains the 
Particle-hole symmetry and the asymmetry slowly becomes obvious with q increasing (the 
first and second panel above in Figjoj). Thus, we can reasonably conclude that whether the 
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FIG. 5: The dispersion relation between k and oj for q = 1.2 (left plot for a = 0.09 and right 
plot for a = -0.19). We fit the data as follows: Co ^ -1.04515(-A;)i-20323 fo^. ^ = o.09 and 
cj ~ — 0.642177(— A;)^''^''^^ for a = —0.19. They show that for q = 1.2, the scaling exponent 5 
decreases with the decrease of the values of GB parameter a and asymptote to 1 for smaller a. 



< 


1 J 




/ 



-3 -1 




FIG. 6: Density plots of the spectral function for T = 0. Lighter regions have higher value of the 
function. The three panels above correspond to a = —0.19 and q = 0, a = —0.19 and q = 0.1, 
a = and q = 1 (from left to right). The three panels below correspond to q = and q = —1, 
a = —0.19 and q = 1, a = —0.19 and q = —1 (from left to right). 
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system has the Particle-hole symmetry or asymmetry is related to the different values of 
q. Some similar conclusions have also been pointed out in the investigations on fermions 
in charged BTZ black hole In addition, we also notice that the Gauss-Bonnet coupling 
constant a has less effects on the (a) symmetry of the system. For comparison, we present 
the plot for a = and g = lorg = — 1 (the last panel above and the first panel below in 
Fig§. 

V. CONCLUSIONS AND DISCUSSION 

We have studied the main features of the fermions in charged Gauss-Bonnet black hole 
for zero temperature limit and massless fermions by AdS / CFT correspondence. The general 
behavior of the spectral function is similar with the case of RN black hole [4J . However, the 
Gauss-Bonnet coupling constant a changes the shape of the spectral function. Especially, 
near the quasi-particle like peak, the effect the a exerts on the shape of the spectral function 
is more significant. Therefore, it is interesting to further understand the behavior of Greens 
functions in the Quasi-topological gravity or the case with Weyl corrections, which has more 
physical contents. 

Furthermore, we especially focus on the dispersion and the Particle-hole asymmetry. 
Their properties can be summarized as follows. Generally, at g = 1, we find a dispersion 
relation uj{k) ~ and the scaling exponent 5 7^ 1, indicating this system is non- Fermi 
liquid. Therefore, the system can be candidates for holographic dual of generalized non- 
Fermi liquids. More importantly, the behavior of this system increasingly similar to that 
of the Landau Fermi liquid when a is approaching its lower bound. It seems that the GB 
term can also describe the dual field theories with both the weak-coupling and the strong- 
coupling. Also we discuss briefly the cases of the charge dependence. At larger values of 
g, new scaling behavior appears. For instance, for certain values of charge q (for example, 
q = 1.2), the scaling exponent 6 decreases with the decrease of the values of GB parameter a 
and asymptote to 1 for smaller a. While for larger charge q (for example, q = 1.5), the scaling 
exponent 6^1 independent of the a. Another important characteristic of this system is 
the Particle-hole asymmetry. From the density plot of the spectral function, we find that 
for zero temperature limit, the Particle-hole (a) symmetry is controlled by the parameter q 
and the Gauss-Bonnet coupling constant a has less effects on it. 
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In addition, other features, such as the log period, studies in Ref. are still preserved, 
which are determined by the geometry AdS2 near the horizon. In the future works, we will 
address the particulars of such properties. In order to searching for the effect of the higher 
curvature corrections on the spectral function of the fermions, we can study the fermions in 
Quasi-topological gravity or the case with Weyl corrections in parallel. 
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Appendix A: Analytical treatment 

In this Appendix, we will give a brief summary on the analytic treatment in the low 
frequency limit developed in Ref.[4lJ- As revealed in Ref.jH], the dispersion relation can be 
given by^ 



In our conventions, gp = 2. The above expressions was derived by the fact that the 
geometry near the horizon is AdS2 x M'^^^. As discussed in SectiorjlH the geometry near the 
horizon is also AdS2 Gauss-Bonnet gravity, independent of a. Therefore, in order to 

obtain the dispersion relation, we only need to work out numerically the Fermi momentum 
kp, which is controlled by UV physics. 

9 The equation (93) in Ref.gl]. 



uj{k) oc k^ , with 




(Al) 



While can be calculated by the equation (57) in Ref. [IT] 




(A2) 
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However, we can also obtain the range of kp by WKB analysis^" 
Ref.|41j. when = 0, the range for allowed kp is^^ 

d-2 _ kp 



< — < 1. 

1) 



. As observed in 



(A3) 



In order to test the robustness of our numerical result on the Fermi momentum. We 
list the values of the Fermi momentum kp determined numerically and the range of kp 

we can see 



determined by Eq.(A3) in Table IV for a = and different q. From Table IV 



that our numerical result is well within the interval determined by Eq. ( A3 ) 



q 


0.5 


1 


1.2 


1.5 


fci? (numerical) 


0.81 


1.8880 


2.349913 


3.0572678 


kp{WKB) 


[0.707, 1.225] 


[1.414,2.450] 


[1.697,2.939] 


[2.121,3.674] 



TABLE IV: The Fermi momentum kp determined numerically and the range oi kp obtained by 
WKB analysis. 
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